Critical Fidelity at the Metal-Insulator Transition 
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Using a Wigner Lorentzian Random Matrix ensemble, we study the fidelity, F{t), of systems at 
the Anderson metal-insulator transition, subject to small perturbations that preserve the criticality. 
We find that there are three decay regimes as perturbation strength increases: the first two are 
associated with a gaussian and an exponential decay respectively and can be described using Linear 
Response Theory. For stronger perturbations F{t) decays algebraically as F{t) ~ t~ 
is the correlation dimension of the Local Density of States. 



where 



PACS numbers: 

The theory of fidelity [l| (also known as Loschmidt 
Echo) has been a subject of intensive research activity 
during the last years (for a recent review see This 
interest has been motivated by various areas of physics, 
ranging from_atomic optics 0,0,01, microwaves [y| and 



^ , microwaves 

elastic waves LZl to quantum information |8|| and quantum 



J3 ' chaos ll [3, 



H, Q El, H [13, [II . It has been 



adopted as a standard measure for quantum reversibility 
and stability of quantum motion with respect to changes 
in an external parameter x^- Formally, the fidelity -F(t), 
is defined as: 



F{t) 



iHct 



1 



(1) 



where Hf and Hb = Hf + a;cB represent the refer- 
ence Hamiltonian and its perturbed variant, respectively, 
while iV'o) is an initial state. One can interpret fidelity 
p]) in two equivalent ways. It can be considered as the 
overlap of an initial state with the state obtained after 
the forward unperturbed evolution, followed by a back- 
ward perturbed evolution. Equivalently, it is the overlap 
of a state obtained after a forward unperturbed evolution 
and the state after a forward perturbed evolution. The 
latter interpretation is closely linked to the concept of de- 
phasing [l9| in mesoscopic devices and coherent manipu- 
lation of a quantum state. Sustaining the coherence of a 
superposition of state vectors is at the heart of quantum 
parallelism in quantum computation schemes [s, 2^, 21 1. 
The first interpretation goes back to the original proposal 
by Peres [H , who used fidelity to study quantum-classical 
correspondence and identify traces of classical (chaotic or 
integrable) dynamics in quantized systems. 

For a quantum system with a classical chaotic counter- 
part, the decay of the fidelity depends on the strength of 
the perturbation parameter x^- Recent studies indicated 
that there are three Xo— regimes: the standard perturba- 
tive regime, the Fermi Golden Rule regime (FGR), and 
the non-perturbative regime. The first two can be de- 
scribed by Linear Response Theory (LRT) leading to a 
decay which depends on the perturbation strength Xr. as 
F[t) ~ e-(^<=*)' and F[t) ~ e"^"*, respectively [13, [ll 
In the non-perturbative regime, the decay is F[t) 



with a rate that is perturbation independent and is given 
by the Lyapunov exponent A of the underlying classical 
system 



Lyapu 



The investigation of the fidelity has recently been ex- 
tended to systems that have integrable classical dynam- 
ics. It was shown [llj that the decay follows a power law 
F[t) ^ t~^'^/'^ , where d is the dimensionality of the sys- 
tem. A similar algebraic decay was found for disordered 
systems with diffractive scattcrers, where now the power 
law is governed by the diffusive dynamics 
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Despite the progress in understanding the fidelity of 
various systems, a significant class was left out of the in- 
vestigation. These are systems which show an Anderson 
metal-insulator transition (MIT) as an external param- 
eter changes. In the metallic regime, the eigenstates of 
these systems are extended, and the statistical proper- 
tics of their spectrum are quite well described by ran- 
dom matrix theory [2I]. In particular, the level spacing 
distribution is very well fitted by the Wigner surmise. 
Deep in the localized regime, the levels become uncorre- 
latcd leading to a Poissonian level spacing distribution 
and the eigcnfunctions are exponentially localized. At 
the MIT, the eigcnfunctions are critical, exhibiting mul- 
tifractal structure characterized by strong fluctuations on 
all scales. The eigenvalue statistics are characterized by 
a third universal distribution [2I, [131 • Representatives 
of this class are disordered systems in d > 2 dimensions, 
two-dimensional systems in strong magnetic fields (quan- 
tum Hall transition), or periodically kicked systems with 
a logarithmic potential singularity [2II . 

Here, for the first time, we address the behavior of F{t) 
for systems at criticality and present consequences of the 
MIT on the fidelity decay. Using the Wigner Lorentzian 
Random Matrix (WLRM) ensemble, we find that there 
are three regimes: (a) the standard perturbative regime 
where the decay is gaussian; (b) the FGR decay where the 
decay is exponential and (c) the non-perturbative regime 
where an initial gaussian decay (Zeno decay) is followed 
by a power law. The latter decay is novel and reflects the 



critical nature of the system. Specifically we found that 

1 



Fit) 



(2) 



where D2 



D2 I d [26[ is the correlation dimension of the 
Local Density of States (LDoS) while is the correla- 
tion dimension of the critical eigenstates and d is the ac- 
tual dimensionality of the system. For the WLRM model 



Dl 



Do since d = 1. The correlation dimension 



D2 is usually defined through the inverse participation 
ratio, P2 ~ J d'^r\il}(r)\'^ ^ L^'^^ , where L is the size 
of the system [23| ■ The correlation dimension is also re- 
lated to the spectral compressibility x = {d—D2)/2d, de- 
fined through the level number variance (SN)'^ w x{^) 



|28l . |29| . |30|. At the same time, D2 manifests itself in 
a variety of other physical observables. As examples 
we mention the conductance distribution 
anomalous spreading of a wave-packet 
dispersion of the diffusion coefficient 34 



31|, I33, the 



the spatial 
3611 , and the 



anomalous scaling of Wigner delay times [37| . 

We use the WLRM model [li,!^!,!!^, defined as: 



H = Ho 



cB 



(3) 



Both Ho and B are real symmetric matrices of size L x 
L with matrix elements randomly drawn from a normal 
distribution with zero mean and a variance depending on 
the distance of the matrix element from the diagonal 



1 



1 



(4) 



where b G (0, L) is a free parameter that controls the crit- 
ical properties of the system (sec Eq. ^ below). Ran- 
dom matrix models with variance given by (ffl were in- 
troduced in [ill and further studied in S lli, [ill, [ill . 
Field-theoretical considerations [H, Isol. l4lj and detail 
numerical investigations [2^ \^ verify that the model 
shows all the key features of the Anderson MIT, including 
multifractality of eigenfunctions and non-trivial spectral 
statistics at the critical point. A theoretical estimation 

si 



for the correlation dimension D2 gives 



/ 46r(3/2)[v^r(i)]- 

1 - 2(27rfo)-i 



1 ,6<1 
,6> 1 



(5) 



where F is the Gamma function. 

The forward and backward Hamiltonians used for the 
calculation of the fidelity (U) are [i^ 



Hf = H(.t) and Hb = H(-.t) 



(6) 



We operate in the basis where Hq is diagonal [40| . In this 
basis, the perturbation matrix B is x~invariant (40| . i.e. 
it preserves the same Lorcntzian power- law shape (j4|. 
while its critical properties (like the multifractal dimen- 
sion D2) remain unchanged. For the numerical evalua- 
tion of F{t), we have used two types of initial conditions 




FIG. 1: Fidelity of an ES, for (a) x = 0.01 (the standard per- 
turbative regime), (b) x = 0.8 (PGR regime) and (c) x = 20 
(non-perturbative regime). The solid lines are the LRT re- 
sults from Eq. (O while the crosses are the outcomes of the 
numerical simulations with model (j3l [ijl . In these simulations 
L = 1000 and b = 10. The mean level spacing of the unper- 
turbed system is set to A « 1. In this case, Xc ~ 0.59 and 
Xprt ~ 1.88. The dotted line in (c) is plotted to guide the eye 
on the power-law behavior. 



iV^o): an cigcnstate of Hq (ES) and a generic "random" 
state (RS). In both cases, the results are qualitatively the 
same. Therefore, we will not distinguish between them. 
In our numerical experiments we used matrices of size 
varying from L = 1000 to L = 5000. We have performed 
an averaging over different initial states and realizations 
of the perturbation matrix B (typically more than 1000). 

An overview of the temporal behavior of the fidelity 
F{t) for three representative perturbation strengths is 
shown in Fig. 1. For perturbation strengths smaller than 

^ ^ A 



1 + T 144[, the decay of F{t) is gaussian (see 



Fig. la). The perturbative border Xc is the perturbation 
strength needed in order to mix levels within a distance 
of a mean level spacing A [i^. Above this border, one 
typically expects an exponential FGR decay of fidelity 
[lO| , with a rate given by the width of the Local Density 
of States (LDoS) ^ (see Fig. lb). We can apply LRT 
[3] to evaluate the decay of F{t) in these two regimes. 
The resulting expression reads 



{Fit))B,no~^-i^xfCit)^C 



(7) 



where (. . ^b.uq represents a double average over B and 
initial states. The right hand side of expression ([7]) as- 
sumes the validity of infinite order perturbation theory. 
The correlator C{t) is 



C{t) 



dri 



dT, 



'C'„(ri-T2)-2Xt2 (8) 
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b = 0.32 
b = 1.00 
b = 3.16 



-t 




the fidelity F(t) clearly displays an inverse power law, 



t 



FIG. 2: (color) Fi for b = 0.32, 1.00 and 3.16. The initial 
state was chosen to be an ES. The mean level spacing of the 
unperturbed Hamiltonian is set to be A « 1, while the size of 
the matrices is i = 5000. In the cases reported here we have 
choose a; = 5. In the inset, we also present the fidelity for RS 
for b — 0.32 and b = 1.00 using the same parameters except 
L = 1000. The straight lines are plotted to guide the eye. 



where T = is the inverse participation ratio of 

the initial state, C„(< - t') = 2(1 + J2j <^ln cos[(£;^*'^ - 

Ei°^){t - t')]), and Ei°^ denotes an eigenvalue of Hq. 
In the case of standard GOE ensembles with ~ 1, 
Eq. ^ reduces to the expression derived in hj. The 
prediction of LRT ([7]) is plotted together with the numer- 
ical results in Fig. [1] for different perturbation strengths. 
A good agreement between Eq. ([7]) and the numeri- 
cal data is observed for perturbation strengths less than 

Ay 



'i.28[7r/2--arctan(i/b)] 



(see Figs. la,b) \M 



i.68[7r/2-arctan(i/6)] 
For X larger than ipit the decay of F(t) cannot be cap- 
tured by LRT (see Fig. Ic). The non-perturbative char- 
acter of this regime was identified already in the frame of 
the parametric evolution of the Local Density of States 
(LDoS) [iOl. A representative temporal behavior of F(t) 
for X > Xprt is reported in Fig. Ic. For short times the de- 
cay of F(t) is gaussian. For longer times, we can observe 
a transition to a power law decay. The initial gaussian 
decay F{t) ~ e~^ *" is universal and can be identified 
with the quantum Zeno effect [l|, Q- K is valid until 
times tz ''^ 1/x. We will focus in the observed power-law 
decay which take place for t > tz- 

The numerical results for three different b- values, b = 
0.32, 1, and 3.16, are reported in Fig. [2l We use the 
time-averaged fidelity 



Fj{t) EE (Fit)), 



F{t)dt 



(9) 



to reduce further statistical fiuctuations. In the inset, we 
present the raw data for the fidelity decay. In all cases 



F{t) cx 



(10) 



with a power 7 that depends on the band-width param- 
eter b. By fitting our data to Eq. (fTU]) . the power-law 
exponent 7 is extracted. In Fig. Owe summarize the ex- 
tracted 7's for both ES and RS initial conditions as a 
function of the bandwidth b. The results are essentially 
identical within the numerical accuracy of our fitting pro- 
cedure. 

If the initial state {ipo) is an eigenstatc of the backward 
(or forward) Hamiltonian then the fidelity is simply the 
survival probability P{t) = \ {ipo\e~^^^'^ \ipo) ^ of wave- 
packet dynamics. In the latter case, it is known that the 
survival probability at criticality decays as P{t) ^ l/t^^ 
[3^. However, in these fidelity experiments, the initial 
state is neither eigenstate of ifb nor of H{. In fact, 
Ref. [l3| shows that the physics of quantum fidelity in- 
volves subtle cross correlations which in general are not 
captured by the survival probability (or the LDoS which 
is its Fourier transform) alone. Motivated by this equiv- 
alence between fidelity and survival probability for the 
specific choice of initial condition \ipo) 1 we have compared 
in Fig. [3] the extracted power law exponents 7 with the 
correlation dimension D2 = -D2 = -^2 The agree- 

ment between the 7 and the D2 is excellent for all b's 
confirming the prediction ([2|). 

The connection between the exponent 7 and the fractal 
dimension calls for an argument for its explanation. 
The following heuristic argument provides some under- 
standing of the power law decay Eq. ([2]). For any fi- 
nite Hilbert space the fidelity F{t) approaches the value 
Foo = 1/i, being the inverse of the dimension of the 
Hilbert space. If the dynamics, however, take place in 
a space with an effective reduced dimension , we will 
have Foo = l/L^? [4^. Assuming a power law decay 
PU)) for the fidelity, we can estimate how the time at 
which F{t^) = Foo scales with L, i.e. ~ L^^/'i, On 
the other hand, the dynamics of a critical system is char- 

acterized by an anomalous diffusive law ^ t^, ^ ^ 
[33! which defines the time ~ L^t/^^ needed to ex- 
plore the available space L. Equating the two expressions 
for we finally get that 7 = D2. Although the numeri- 
cal results leaves no doubt on the validity of Eq. ^ , a 
rigorous mathematical proof is more than desirable. 

In conclusion, we have investigated the fidelity de- 
cay for systems at MIT. Depending on the perturbation 
strength x, we have indentified three distinct regimes: 
For x < Xc the fidelity decay is gaussian; for Xc < x < 
Xpj-t the decay is exponential and for x > Xp^t the decay is 
power law. The first two regimes are described by LRT. 
The third is non-perturbative. The power law decay is 
dictated by the critical nature of the system. Specifically, 
we have found that the power-law exponent is equal to 
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FIG. 3: (color) The fitting parameter 7 for b — 
0.03, 0.10, 0.32, 1.00, 3.16, and 10. We are using A = 1, x = 5. 
The analytical (solid lines) Eq. ((5)1 and numerical (crosses) 
[45| results for D2 are also shown for comparison. 



the correlation dimension of the critical eigenstates. 
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